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SYNOPSIS 

This  report  presents  two  methods,  consistent  deformations  and  an  energy 
method,  for  the  stress  analysis  of  a  girder-arch  bridge.   Using  the  first 
method,  two  solutions  were  calculated,  the  first  of  which  takes  into  account 
the  axial  stress  effects  in  the  structure,  and  the  second  ignores  these 
effects.  A  third  solution  is  developed  using  the  energy  method. 

For  the  method  of  consistent  deformations,  the  first  step  is  to  find  the 
deflections  of  the  girder  and  of  the  arch  separately.   The  redundant  forces  in 
the  vertical  member  are  then  found  by  requiring  the  deformations  of  the  girder, 
the  arch  and  the  vertical  members  to  be  consistent.   The  moment  influence 
lines  can  then  be  drawn.   The  second  method,  which  uses  a  minimum  energy 
approach,  treats  the  structure  as  a  single  system  by  assuming  that  the  por- 
tion between  the  girder  and  the  arch  acts  as  one  member.   Using  the  virtual 
displacement  method,  the  redundant  reactions  of  the  structure  can  then  be 
found,   fiext,  the  total  moment  at  each  section  of  the  compound  structure  can 
be  calculated.   Then  the  bending  moment  in  the  girder  and  In  the  arch  can  be 
found  for  each  section  by  distributing  the  total  moment  induced  in  the 
compound  section  between  the  two  elements.   Using  these  results  the  moment 
influence  lines  for  the  girder  and  the  arch  can  be  drawn. 

The  computations  in  the  report  are  generally  expressed  in  matrix  form 
and  have  been  wholly  performed  by  a  computer  so  that  the  tediousness  and 
possible  inaccuracy'of  hand  computation  are  minimized.  A  comparison  of  the 
three  solutions  Indicates  that  the  results  are  reasonably  close  to  each  other. 


INTRODUCTION 


.1 


The  girder-arch  bridge  (Flg.l)  analysed  In  this  report  is  a  mordern  type 
of  bridge  which  combines  the  girder  and  the  arch  by  using  vertical  hinged 
members  to  connect  the  two  main  members.   Sy  combining  the  girder  and  the  arch, 
the  stiffness  of  the  system  is  greater  than  that  of  a  simple  girder  or  a 
simple  arch.   The  moment  induced  in  the  system  is  greatly  reduced  by  the 
action  of  the  arch.   Furthermore,  girder-arch  structures  are  economical  (1)" 
and  are  pleasing  In  appearance.   This  type  of  bridge  is  especially  suitable 
for  use  as  a  long  span  over  a  deep  valley. 

Of  the  many  possible  methods  of  analysis  for  this  type  of  bridge,  this 
report  presents  two  which  may  be  used  to  analyse  high-order  statically  Inde- 
terminate structures.  A  digital  computer  was  available  and  therefore  the 
calculations  for  both  of  these  methods  were  wholly  performed  by  the  digital 
computer. 


Fig.  1.   Elevation  of  a  girder-arch  bridge. 


Numerals  in  parentheses  refer  to  corresponding  Items  In  the  References. 


Since  the  bridge  analysed  in  this  report  has  variable  cross  sections  in  the 
arch,  this  report  shows  a  general  procedure  for  a  computer  analysis  to  handle 
a  structure  with  variable  cross  sections.   Because  computers  are  now  readily 
available,  either  of  the  methods  presented  in  this  report  can  be  used  to 
handle  most  structures  with  variable  sections  such  as  continuous  beams  or 
rigid  frames  very  easily.  Thus  the  methods  this  report  presents  are  not 
restricted  to  this  particular  problem  but  can  be  generalized  to  analyse  many 
indeterminate  structures. 


GENERAL  ANALYSIS  PROCEDURES 

I.  Consistent  Deformation  Method. 
a.  Solution  1. 

Considering  both  flexural  and  axial  stress  effects. 

This  method  requires  that  the  deflections  of  the  girder,  the  vertical 
truss  members,  and  the  arch  be  consistent  at  each  node  point  when  the  structure 
is  loaded.  Therefore  the  first  step  is  to  find  the  deflections  of  the  girder 
for  each  node.   For  drawing  the  influence  lines,  the  load  acting  along  the 
structure  is  taken  as  a  moving  unit  load.   This  step  is  accomplished  by  using 
the  conjugate  beam  method.   The  next  step  is  to  compute  the  deflections  of  the 
two-hinged  arch.   For  this  purpose,  the  deflections  of  the  simple  arch  are 
first  calculated,  then  the  superposition  method  is  applied  to  find  the  actual 
horizontal  reactions  of  a  two-hinged  arch.   Knowing  the  horizontal  reactions 
and  using  the  method  of  superposition  again,  the  vertical  deflections  of  a 
two-hinged  arch  are  obtained.   The  third  step  is  to  set  up  the  equations  for 
calculating  the  redundant  forces  in  the  vertical  members.   Knowing  the  deflec- 
tions of  the  girder  and  of  the  arch  due  to  unit  load  acting  along  the  girder 
and  the  arch  respectively  from  the  first  and  second  step,  the  consistent  de- 
formation relations  can  then  be  applied  for  each  node  point.   Since  the  con- 
sistent deformation  relationships  require  that  the  deflection  of  the  girder 
at  each  node  point  equal  to  the  deformation  of  the  vertical  member  plus  the 
deflection  of  the  trwo-hlnged  arch  at  the  same  node  point,  these  relationships 
can  be  set  up  for  each  node  point  as  a  function  of  the  vertical  member  redun- 
dant forces.   The  vertical  redundant  forces  can  then  be  obtained  by  solving 
these  simultaneous  equations.  With  the  redundant  vertical  member  forces  known, 


it  is  possible  to  calculate  the  redundant  horizontal  reactions  or  the  whole 
structure  by  superposition.   Knowing  all  the  redundant  forces  acting  on  the 
structure,  the  bending  moment  for  each  section  of  the  girder  and  the  arch 
can  be  calculated.   The  steps  required  in  this  analysis  can  be  outlined  as 
follows: 

(1)   Computation  of  the  deflections  of  the  simple  girder. 
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Fig.  2A.   Notation  for  node  points  of  a  simple  girder. 

Subdividing  the  girder  with  node  points  and  placing  the  loading,  P,  =  1, 
at  node  point  j  as  shown  in  Fig.  2A,  the  deflections  of  the  beam  at  any 
section  1<  can  be  found  by  using  the  conjugate  beam  method. 

For  the  loading  shown  in  Fig.  2A,  the  M/EI  diagram  is  shown  in  Fig.  2B. 

PjjX(m  -  j)X 


Fig.  2B.  Loaded  conjugate  beam. 


This  M/EI  diagram  is  the  loading  diagram  for  the  conjugate  beam  as  shown  in 
Fig.  2B.  The  left  hand  reaction,  R,,  of  this  conjugate  beam  can  be  found  by 
taking  moments  about  B: 

.       F,(jX)<m  -  j)X   .        ,. 

K    =  —  J ±r  fir  +  <ra  -  Utf 
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Since  the  moment  at  a  section  in  the  conjugate  beam  is  the  deflection  of  the 
real  beam  at  that  section,  the  deflection  at  any  section  k  of  the  real  beam 
can  be  found  by  taking  moments  about  node  point  k  of  the  conjugate  beam. 
For  k  <  j : 
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For  k  >  j  : 


.    _  L3j(m  -  k)[k(2m  -  k)  -  js1 

6:<j  =  TC (1"b) 

on   n.1  • 


The  first  and  second  subscripts  denote  section  and  loading  point  respectively 
and  are  used  throughout  this  report  unless  otherwise  specified. 


(2)   Computation  of  the  deflections  of  the  two-hinged  arch. 
Using  the  method  of  virtual  work,  the  vertical  deflection  of  any  section 
k,  with  loading,  P.  =  1,  acting  at  node  point  j,  as  shown  in  Fig.  3  can  be 


ds  = 


dx 
cosS 


Fig.  3.  Moment  diagram  of  a  simple  parabolic  arch  with 
variable  cross  sections. 


found  by  the  following  formula  (2) (3). 


L  m,m,dx      L  n.n^dx 
sr     .   J  »<       .    J  « 

6ki  =  i  ~ —  +  s  ~ 

J     0  EI  cose     0  A  Ecos6 


(2) 


If  the  load  Ffc  is  acting  horizontally  at  hinge  point  D,  the  corresponding 
values  of  moment  m,  and  thrust  n.  are  y  and  cos0  respectively.   Thus  if  y  and 
cos9  are  substituted  for  m.    and  n.  respectively  In  equation  (2),  the  horizon- 
tal displacement  of  point  D  due  to  the  vertical  load,  P.  =  1,  can  be  found. 
Furthermore,  if  y  and  cose  are  substituted  for  m,  and  n  respectively  In 
equation  (2),  the  value  of  6.  .  is  the  horizontal  displacement  of  point  D  due 
to  a  unit  horizontal  load  acting  at  point  D.   If  the  thrust  effects  are  neg- 
lected, equation  (2)  can  be  reduced  to 
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Knowing  the  deflections  of  a  simple  arch,  the  superposition  method  (Fig.  4) 
can  now  be  used  to  find  the  deflections  of  a  two-hinged  arch. 


Fig.  4.  Deflections  of  a  two-hinged  arch. 


With  the  loading,  P.  =  1,  acting  at  node  point  j,  the  expression  for 
the  deflections  at  node  point  k  can  be  expressed  as  follows  (3): 
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It  can  be  seen  that  the  horizontal  reaction  H  in  equation  (5)  also 
comes  from  equation  (4)  by  setting  the  left  hand  side  of  the  equation  equal 


to  zero.  This  is  simply  the  superposition  method  for  calculating  the  hori- 
zontal reaction  of  a  two-hinged  arch.   The  numerator  in  the  right  hand  side 
of  equation  (5)  represents  the  horizontal  displacement  due  to  a  unit  vertical 
load,  P.  =  1,  acting  at  point  j,  and  the  denominator  represents  the  horizontal 
displacement  due  to  a  unit  load  acting  horizontally  at  point  D.   These  inter- 
pretations can  be  made  based  on  equation  (2)  as  stated  before. 

If  the  influence  of  thrust  in  the  arch  is  neglected,  equation  (5)  can 
be  reduced  to 


(6) 

o    Tcose    ' 

(3)   Calculation  of  the  redundant  forces  In  the  vertical  members. 

Fig.  5  shows  the  whole  system  of  a  girder-arch  bridge.  For  the  loading, 
P.  =  1,  acting  at  j,  there  are  forces  induced  in  every  vertical  member  as 
shown.   From  the  relations  of  consistent  deformations,  the  deflection  of  the 
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Fig.  5.   Redundant  forces  in  the  ties  of  a  girder- 
arch  bridge. 

girder  at  each  node  point  must  be  equal  to  the  deformation  of  the  vertical 
member  plus  the  deflection  of  the  arch  at  the  same  node  point.   This 
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relationship  can  be  expressed  as: 


6k  =  «k  +  6i<  •  (7> 


a 

The  deflections  of  the  simple  girder,  6",  and  of  the  two-hinged  arch, 

r 
6.  ,  have  already  been  calculated  by  the  use  of  equations  (1)  and  (4)  res- 

v 

pectively,  and  the  deformation  of  the  vertical  member,  6. ,  is  given  by: 

v  _  XkLk 

k     "     AVE  • 


With  the  loads  acting  on  the  structure  as  shown  in  Fig.  5,  it  may  seen  that 
the  various  deflection  at  point  1  can  be  expressed  as: 


61  =  5I,1(P1  "  V  +  61,2(P2  "  V  +  -  +  6l,n(Pn  "  V  > 


v     X1L1 

6i  =  77, 


h  -  °  • 


Substituting  these  values    into  equation   (7),    it  follows   that: 

ll(Pl   "  V   +  6?,2<P2    "  V   +  -  +  6l,n<Pn  "  V    "  ^  M 

A  E  • 


Similarly,  the  deflections  at  point  2  can  be  expressed  as: 


62  -  62,1(P1  "  V  +  62,2(P2  "  V  +  •"  +  62,n(Pn  "  V  • 

v     X2L2 

62  =  -v- 
1  A  E  , 
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62      "     °  +  62,2X2   +  62,3X3   +  -    +  62,n-lVl   +  ° 


Substituting  these  values    Into  equation   (7),    it  follows   that: 


6°     (P,    -  X)   +  6°      (P,   -  Xj   +  •••   +  6,      (P     -  X  ) 
2,1     11  2,2     2         2  2,n     n         n 


+  5,   ,X     +  6,     X     +  •••   +6  X  (b) 

.v„  2,2   2  2,3   3  2,n-l  n-1 


Similar  expressions  can  be  derived   to  express   the  deflection  relation- 
ships  for  the  other  node   points.      Rearranging  equations    (a),(b),...,    it 
follows   that: 

(5l,l+^)X1   +  5f,2X2   +  -+6I,nXn 

"     6I,1P1  +  6!,2P2  +  -   +6l,nPn,  (a)" 


and 


L2 
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2,1    1  2,2  2,2        ,VE     2  2, n-1  2, n-1     n-1  2,n  n 

-     6*   ,P,    +  6?     P     +  •••   +  bl     P  (b)' 

2,11         2,2   2  2,n  n  , 


Expressing  equation's    (a)*,    (b)' ,    ...    in  matrix  form,    it   follows  that: 


12 


- 

M  _ 


+ 

c 

60    C 


+ 


^ 


Ed 

> 


+ 
OJ 

bo  oj 


%c 


+ 


°?oc 


13 


c 


C 


C 
60    C 


14 


or  symbolically 

MM  ■  [BIp]  • 


and 


[X]  =  [A]     [B][P] 


(9) 


(10) 


For  finding  the  influence  lines  for  the  girder  and  the  arch,  the  [P"l 
matrix  is  taken  as  a  unit  matrix. 

(4)   Calculation  of  redundant  horizontal  reactions. 

It  may  be  seen  from  Fig.  4  that  whenever  a  load  acts  at  j,  there  is  a 
corresponding  horizontal  reaction  H.,  which  has  been  found  from  equation 
(5).  From  Fig.  5,  it  may  be  also  seen  that  whenever  a  system  of  loads  acts 
on  the  structure,  there  is  a  corresponding  system  of  vertical  member  forces 
X, ,  X,,  ...,  X   induced  in  the  vertical  members.   The  horizontal  reactions 
corresponding  to  such  "sots"  of  vertical  redundant  forces  as  calculated  in 
the  preceding  section  can  be  expressed  as: 


h!  >  x,.h!|  +  x  h*  + 

J      2j  2      3]  3 


+  X   ,  ,H  .   . 
n-l,j  n-1 


(U-a) 


With  the  unit  loading  acting  along  the  girder,  equation  (11-a)  can  be  ex- 


pressed in  matrix  form  as  follows: 
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(11-b) 
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or  symbolically 

[/]  =  [X'jH6]  .  (11-c) 

It  may  be  seen  that  the  matrix  [X*  "1  Is  the  transpose  of  matrix  [X^  obtained 

from  the  preceding  section  if  the  first  and  the  last  rows  are  eliminated. 

(5)  Calculation  of  moment  influence  lines  for  the  girder. 
Knowing  the  redundant  forces  induced  in  the  vertical  members,  it  is 
readily  seen  from  Fig.  5  that  the  moments  in  the  girder  can  be  expressed  as: 
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If  point  1  is  allowed  to  vary,  equation  (12-a)  can  be  expressed  in  matrix 
form  as: 
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or  symbolically 


pH     =     [mg][Qs]     . 


(12-c) 


(6)   Calculation  of  moment  influence  lines  for  the  arch. 

Using  a  line  of  reasoning  similar  to  that  explained  in  the  preceding 
section,  it  can  be  seen  that  the  final  moment  influence  line  for  the  arch 
can  be  expressed  as : 
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or  symbolically 

[M1"]  =   [»r][x] 


(13-b) 


b.  Solution  2. 

Considering  the  flexural  stress  effects  only. 

This  solution,  which  uses  the  concept  of  consistent  deformations  but 
neglects  the  axial  stress  effects,  is  quite  similar  to  solution  1  which  con- 
sidered both  the  flexural  and  axial  stress  effects.  With  the  axial  stress 
in  the  vertical  members  neglected,  it  may  be  seen  that  the  girder  can  not 
deflect  at  points  1  and  n,  therefore  the  girder  is  now  considered  as  a  con- 
tinuous girder  with  supports  at  A,  1,  n  and  B  as  shown  in  Fig.  6. 

To  apply  the  consistent  deformations  relation  at  each  node  point,  it 
is  first  necessary  to  calculate  the  deflections  of  the  continuous  girder  for 
each  node  point.   To  calculate  the  deflections  of  the  continuous  girder  for 
each  node  point,  it  is  necessary  to  know  the  redundant  reactions  at  supports 
1  and  n.   This  is  accomplished  by  using  the  superposition  method,  taking 
advantage  of  the  known  values  of  the  deflections  of  a  simple  girder  as  com- 
puted in  solution  1.   Knowing  the  values  of  the  redundant  reactions,  the 
superposition  method  is  used  again  to  obtain  the  deflections  of  the  con- 
tinuous girder. 

The  second  step  is  to  find  the  deflections  of  a  two-hinged  arch.   Since 
the  axial  stress  effects  are  neglected,  equations  (3)  and  (6)  are  used 
instead  of  equation  (2)  and  (5).  Substituting  the  values  from  these  two 
equations  Into  equation  (4),  the  deflections  of  a  two-hinged  arch  can  be 
calculated  as  stated  previously  in  solution  1.   Knowing  the  deflections  of 
the  continuous  girder  and  of  the  two-hinged  arch  for  each  mode  point,  the 
relationship  of  consistent  deformations  can  be  applied  to  each  node  point, 
and  the  equations  for  the  unknown  vortical  member  forces  can  be  set  up. 
Solving  the  equations  simultaneously  and  calculating  the  redundant  horizontal 
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reactions  as  done  in  solution  1,  the  redundants  in  the  structure  can  be 
completely  solved  and  the  moment  influence  lines  can  be  drawn.   The  steps 
involved  in  this  analysis  can  be  outlined  as  follows: 

(1)  Calculation  of  the  deflections  of  the  continuous  girder. 

Since  the  axial  stress  effects  are  neglected,  the  girder  is  considered 
as  a  continuous  girder  with  supports  at  A,  1,  n  and  B  as  shown  in  Fig.  6. 


Fig.  6.  Superposition  of  the  deflections  of  the 

continuous  girder. 


Using  the  deflections  of  the  simple  girder  as  calculated  from  equations 
(1-a)  and  (1-b)  in  solution  1,  the  deflections  of  the  continuous  girder  can 
be  found  by  using  the  superposition  method  as  shown  in  Fig.  6.   Since  the 
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deflections  of  the  girder  at  points  1  and  n  have  to  be  zero,  it  follows  that: 


8,  ,R,  .  *  6,  R  .  -  6,  .  =0  , 
"1,1  I,:  l,n  n,j    l,j 


6  ,R,  .  +  6   R  ,  -  8  :  *  0  . 

n,l  l,j    n,n  n,j    n,j 

Solving  these  two  equations  simultaneously,  it  follows  that: 


6,  .8    -  6   .6, 
l,j  n,n    n,j  l,n 
R,  .  = 


l'J    5,  ,6    -  6,   6  , 

1,1  n,n    l,n  n, 1   , 


(14) 


R^  • 


n,j  1,1    l,j  n,l 


n'J    6,  ,6    -  6,   8  , 
1,1  n,n    l,n  n, 1 


Knowing  the  values  of  R,    .    and  R      ,,    the  deflections  of   the   continuous  girder 
l»J  n,j 


6      .     =     6Svg.   -  R,    .6fS,    -  R      .8»*  (15) 

k,j  k,j  l,j    k,l  n,j    k,n      . 


(2)  Calculation  of  the  deflection  of  the  two-hinged  arch. 

This  step  is  exactly  the  same  as  stated  in  step  (2)  of  solution  1, 
Equations  (3),  (4),  and  (6)  are  used  in  this  solution  whereas  equations  (2), 
(4),  and  (5)  are  used  in  solution  1. 

(3)  Calculation  of  the  redundant  forces  in  the  vertical  members. 

This  step  can  be  seen  to  be  similar  to  step  (3)  in  solution  1.  With 

v 
the  axial  stress  effects  ignored,  the  term  6^  in  equation  (7)  is  omitted. 

Noting  that  in  this  solution  there  are  no  deflections  at  points  1  and  n,  it 

may  bo  seen  that  the  first  and  the  last  rows  of  the  first  matrix  in  the  left 

hand  side  of  equation  (8)  are  omitted.   Furthermore,  noting  that  the  redun- 

dants  X.  and  X  which  are  the  values  of  R,  .  and  R   .  have  been  calculated, 

1      n  l,j      n,j 
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equation   (S)   can  be  expressed,    in  this   solution,   as: 


62,2+62,2 
63,2+63,2 


6!,3+62,3 


6,   ,+6,  , 
o,3     3,3 


62,n-2+62,n-2 


63,n-2+53,n-2 


*g 


"-s: 


2, n-1     2, n-1 


63,n-l+63,n-l 


6n-2,2+6n-2,2      6n-2,3+5n-2,3 


6s         +8'  8s         +6 

n-1,2     n-1,2       n-1,3     n-1,3 


62,2  52,3 


63,2  63,3 


fig  Kg 

n-2, 2     °n-2,3 

6n-l,2       n-1,3 


.+6 


t» 


+6r 


n-2,n-2     n-2,n-2        n-2,n-l     n-2,n-l 


6s  +6  6s  +8 

n-l,n-2     n-l,n-2       n-l,n-l     n-l,n-l 


tS  *§ 


n-2 


n-1 


2, n-2  2, n-1 


3, n-2  3, n-1 


5n-2,n-2      6n-2,n-l 


n-1, n-2        n-1, n-1 


n-2 


n-1 


(16) 


or  expressed   symbolically   in  the   same   form  as  equation   (9): 
[A][X]    =     [b][p1    , 


and 


[X]     =      [A]     [B][P1     . 


(17) 


(18) 


(4)  Calculation  of  the  redundant  horizontal  reactions. 

With  a  line  of  reasoning  similar  to  that  explained  in  step  (4)  of 
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solution  1,  and  noting  that  whenever  the  loads  act  at  points  1  and  n  they 
will  cause  no  stress  in  the  remaining  vertical  truss  members,  It  may  be  seen 
that  the  first  and  the  last  rows  in  the  matrices  [H  1  and  [X*  ]  of  equation 
(11-c)  are  eliminated.   Thus  equation  (11-b)  can  be  expressed,  in  this 
solution,  as: 


n-2 
n-1 


r 


"2,2 


3,2 


2,3     3,3 


X2,n-2  X3,n-2 


X2,n-1  X3,n-1 


n-2, 2     n-1, 2 


X 


n-2, 3     n-1, 3 


X 


n-2, n-2  '  n-1, n-2 


X        X 

'  n-2, n-1   n-1, n-1 


t 

H2 

4 

• 

• 

• 

Hn-2 

[«n-l_ 

(19) 


or  symbolically, 


[Hfl  =  [k'jH*] 


(20) 


It  may  be  seen  that  the  matrix  fXf  ~|  is  the  transpose  of  the  matrix  TX"] 
obtained  from  the  preceding  section.   Therefore  equation  (20)  can  be  ex- 
pressed as: 

T 


[Hfl  =  pcfpr*] 


(21) 


(5)  Calculation  of  moment  influence  lines  for  the  girder. 

In  order  to  draw  the  moment  influence  lines  for  the  girder,  the  in- 
fluence values  of  the  moments  induced  in  the  continuous  girder  for  each 
section  due  to  a  unit  load  acting  along  the  girder  must  be  known.  With  the 
values  R    and  E   .  from  equation  (14),  the  reactions  at  the  end  points 
can  be  found  by  taking  moments  with  respect  to  either  end  of  the  girder 
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(Fig.    6)   and  can  be  expressed  as: 

RA,j  "      [(ra  "   *>*  "  Rlj(L2  +  L3}    "  Rn,jL2>L     • 


(22) 


Knowing  the  values  of  R.   and  R   ,  the  moments  at  each  section  with  the 
loading,  P.  =  1,  acting  at  j  can  be  expressed  as: 


for  i  <  j 


m    =  R.  .(IX)  +  R,  .(i  -  1)X 


for  1  >  J   : 


(23-a) 


m.  .   =  RE  .(in  -  i)X  +  R   .  (n  -  i)X 

i.J      B,j  n,j 


(23-b) 


Knowing  all  the  values  of  m.  . ,  the  same  arguments  used  in  solution  1  for 

i»  J 

finding  the  moment  influence  lines  can  be  applied  in  this  solution.  Noting 

that  when  the  loads  act  at  points  1  and  n  they  will  cause  no  moment  at  any 

g 
section  of  the  structure,  the  first  and  the  last  columns  of  the  matrix  [M  "1 

and  m   in  equation  (12-c)  are  eliminated.   Equation  (12-b)  can  then  be  ex- 
pressed, in  this  solution,  as: 


MS      M8 
1,2    ™1,3 


M2,2    %3 


"g  1  ■>     "*  i  ■» 
n-1,2   n-1,3 


M8  ,    MS  , 
n,2     n,3 


l,n-2     l,n-l 


2,n-2     2,n-l 


n-l,n-2   n-l,n-l 


MS   ,    MS   . 
n,n-2     n,n-l 


24 


1,2 

4. 


■!,3 

m2,3 


n-l,2        n-1,3 


n,2  n,3 


"U-2         "In-l 


"2,n-2         '"2,n-l 


i8   ,        .,     mS   , 
n-l,n-2        n-l,n-l 


n,n-2  n,n-l 


UX2,2        "X2,3   •"    'X2,n-2        "  X2,n-1 


X3,2        1_X3,3    "'    "X3,n-2        "  X3,n-1 


n-2,2  "     n-2,3* 


Vl,2   "  Xn-1,3' 


•  1-X     _       _   •  X     -        . 
n-2,n-2         n-2,n-l 


X      ,        „    1-X      , 

n-l,n-2         n-l,n-l 


(24) 


or  symbolically. 


[MS1     =     [n-g][Q?T 


(25) 


(6)  Calculation  of  moment  Influence  lines  for  the  arch. 

With  the  same  interpretations  as  those  stated  in  step  (6)  of  solution  1, 
and  noting  that  there  are  no  deflections  or  moments  when  the  loading  acts  at 
points  1  and  n,  equation  (13-a)  can  be  expressed,  in  this  solution,  as: 


•£,2    ^,3 


M3,2    M3,3 


2,n-2     2,n-l 


M,   -    M,    , 
3,n-2     3,n-l 


M  ,  ,  M  ,  ,  •  •  •  M  ,   ,  K   ,    , 
n-2,2   n-2,3        n-2,n-2   n-2,n-l 


Mr  .  ,  M37  ,  ,  •  •  •  M1"  ,    ,  Mr  ,    , 
n-1,2   n-1,3        n-l,n-2   n-l,n-l 
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2,2  2,3 

"3,2  m3,3 


2,n-2  2,n-l 


3,n-2  3,n-l 


mn-2,2     ran-2,3   '    *   *  V2,n.2     mn-2,n-l 


r  r  r  r 

n-1,2        n-1,3  n-l,n-2        n-l,n-l 


X  X  •    •    •  X  X 

2,2  2,3  2,n-2        2,n-l 


X3,2  X3,3    '    "    '   X3,n-2     X3,n-1 


Xn-2,2     Xn-2,3   "'  Xn-2,n-2Xn-2,n-l 


Xn-1,2     Xn-1,3    ""    Xn-l,n-2Xn-l,n-l 


(26) 


or  symbolically, 


(27) 
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II.  Energy  method. 

Since  the  structure  analysed  in  this  report  is  14  degrees  statically 
indeterminate,  the  ordinary  energy  method  which  requires  writing  the  energy 
terms  for  each  segment  of  the  structure  seems  too  tedious.   The  problem  is 
treated  here  by  using  a  method  suggested  by  Professor  James  Michalos  (4). 
Using  this  method,  the  portion  l,C,D,n  of  the  structure  shown  in  Fig.  7  is 
considered  as  a  compound  section  with  flexural  rigidity,  EI,  equal  to  the 
sum  of  that  of  the  girder  and  of  the  arch  in  the  X-axis  direction.   To  find 
the  moments  in  the  girder  and  in  the  arch,  it  is  necessary  to  find  the  reac- 
tion redundants  and  to  calculate  the  moments  induced  in  each  compound  section. 
Then  using  equation  (39),  which  will  be  discussed  later,  the  total  moment  is 
distributed  between  the  girder  and  the  arch.   The  moment  Influence  lines  can 
then  be  drawn.   The  steps  required  in  this  analysis  can  be  outlined  as 
follows: 


C  D 

Fig.  7.  The  reduced  structure  of  a  girder-arch  bridge. 


(1)  Calculation  of  the  redundant  reaction  Influence  lines  for  the 
girder-arch  bridge. 
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The  method  of  virtual  displacements  (5)  Is  used  to  find  the  Influence 

lines  for  the  reactions  V  ,  V„  and  H  as  shown  in  Fig.  7.   Since  the  theory  of 
c   D 

virtual  displacements  states  that  the  deflection  curve  of  a  structure  due  to 
unit  displacement  In  the  direction  of  the  unknown  force  represents  the 
Influence  line  of  that  unknown  force,  the  influence  lines  for  V  ,  V  and  H  can 
be  found  as  follows: 

The  moments  caused  by  the  three  redundant  reactions  V  ,  V_  and  H  are 
shown  in  Fig.  8  (b),  (c),  and  (d)  respectively.   The  total  moment  Induced  in 
the  structure  by  these  redundant  reactions  Is  the  sum  of  the  three  moments. 
Knowing  the  total  moment,  which  is  a  function  of  V  ,  V_  and  H,  for  each 
section,  it  may  be  seen  that  the  total  internal  strain  energy  U  can  be 
expressed  as: 

M  (!"c  +md)Sdx    L!+L2  (mc  +md  +mh^dx    L     <mc  +rad)2dx 
0      2EIS        Lj      2E(Ig  +  ircos9)    Lj  +L2    2EIS 

I,  (Vcmx   +  VDmn)3dx    Ll+I,2  (V^  +  V0mn  +  Hy)2dx 

=     /       +  /  

0  2EIg  Lj  2E(Ig  +  ircose) 

+  /    : (28) 

Lt+L2       2EIb 

According  to  Castlngliano's  first  theorem  (2),  the  first  partial  derivative 
of  the  total  strain  energy  of  the  structure  with  respect  to  the  applied 
action  denotes  the  deflection  component  of  the  point  of  application  of  an 
action  on  a  structure,  in  the  direction  of  that  action.  Or  expressed  mathe- 
matically: 
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(c) 


Moment  diagram  due  to  reaction  V 


Moment  diagram  due  to  reaction  V 


(d)  Moment  diagram  due  to  reaction  II . 


Fig.  8.  Moment  diagrams  due  to  external  loading  and 
redundant  reactions. 


= 

Ac 

9 

-au 

= 

AD 

i 

an 

= 

AH 

■ 

(29) 
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Substituting  equation  (28)  into  equation  (29),  it  follows  that: 


a     .  *H  .  -£r/  y°idx    -,  +  Vr  i+L2       yy1* 

H     "  0H    "   E    4        (r«  +  ircose)    +   E   t,         (ig  +  ircose) 


H   rfLl+L2  ysdx  -,     , 

+  e[/      i — ; — ]  (30-a) 

Lj  (Ig  +  loose) 


T>U            vr      Ll  (mi>Sdx          4+L,           (m.)2dx                 L          (m.)sdx 

A       =  —  =  —  [7     — +  I +  /        — "1 

C            *VC           E     0  j8            L,             (i8  +  ircose)        L1+L2           J* 

+   !°  [/Ll  "I"'"''*  +  ;Ll+L2          mlmndx            +    ,L          mlmndx-1 

E     0  l«            Lj             (!«  +  iroosP)        L,+L2        Is 


„       Li+L2  ym.dx 

Ll  (I8  +  IcosP) 


D  DVn    "      E     0  jg        +Lj  (I«  +  ircose)    +  L,+L2         |« 

v_       L.    (m   )=dx  L.+L,  (m   )sdx  L  (m   )Edx 

♦TC     -^T-*/         5 +  '        — 1 

E     0  18  Lt  (1*5  +  Ircose)        L^  Is 

+  ~  [/  E 1      (30-c) 

L     4  (i8  +   jrcosfi)      f 

or   symbolically, 
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EA       =     E 


I  Vc  +  b  VD  +  c  H 


•2U 


IV     +eV„   +  f  H 
C  D 


(31) 


EA        =  E 

D  -zv 


g  Vc   +  h  VQ  +  1  H 


or   in  matrix  form, 


AH~ 

—           - 

-a  U 
-a  h 

Ac 

=     E 

-au 

= 

AD 

•z>u 

u        — 

-      [A] 


(32) 


Using  the  theory  of  virtual  displacements,  the  Influence  values  of  V  , 
V_  and  H,  which  are  used  to  draw  the  influence  lines,  can  be  found  by  setting 
A„,A„  and  A„,  In  turn,  equal  to  unity  while  the  remaining  terms  in  the 
same  set  of  displacements  is  held  equal  to  zero,  I.e. 


0     0 


[A] 


(33) 


The  subscripts  1,2  and  3  represent  the  order  of  sets  of  the  values  of  V,  Vn 
and  H,  the  total  moment,  Mkt.,   can  be  expressed  as: 


let 


vc  mki   +  vn  "V-n   +  H   ?!, 


1=1.2,3. 


(34) 
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The  deflections  of  the  structure  can  then  be  found  by  using  the  method  of 
virtual  work  and  can  be  expressed  as: 

i     rLi  MktVx    h+L2     MktVx    ,.  ,L    Mktmkdx 

6k   =  *         Z +  '  1 r J  E 

0     EIS     Lt     e(I°  +  I  cose)   Ll+L2   EI      ' 

1=1,2,3.  (35) 

Equation  (35)  gives  three  sets  of  values  of  deflections  which  are  the  in- 
fluence lines  for  H,  Vc  and  V  respectively. 

(2)  Calculation  of  moment  influence  lines  for  the  girder  and  the  arch. 

Knowing  the  influence  lines  for  V  ,  V  and  H,  the  total  moment  induced 
In  each  section  of  the  structure  when  the  structure  is  loaded  can  be  found 
by  adding  the  four  moment  diagrams  as  shown  in  Fig.  8.   The  total  moment,  M  , 
can  be  distributed  between  the  girder  and  the  arch  by  the  following  method  (4). 

The  structure  will  deform  when  the  loads  act  on  it.   If  the  angle  changes 
along  the  girder  and  the  arch  are  equal  at  each  corresponding  point,  that  Is, 
If: 

(36) 


EI  cos0      EIa   , 

where  6  is  the  angle  between  the  X-axis  and  tangent  of  the  rib  at  each  node 
point.  Then 

Mr     Ireose 


(37) 


If  we  denote  M  as  the  total  moment  that  is  resisted  by  bending  in  the  arch 
rib  and  girder  at  the  section,  that  is, 

M   =  Mg  +  Mr   ,  (38) 
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then  from  equation  (37), 


Ms  =  M  ( 1 )   , 

C  Tg 


Is  +  I  cose 

(39) 

i^cose 


if      =  H<      IcosB   )   . 


t   2 

I5  +  I  COS? 


Where  M  can  be  found  by  the  following  equation  as  stated  previously: 

Mt  "  mlj  "  Vil  -  Vin  "  Hyi  .  (40) 

Equation  (39)  thus  provides  a  simple  method  for  drawing  the  moment  influence 
lines  for  the  girder  and  the  arch. 
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NUMERICAL  EXAMPLE 

All  of  the  calculations  for  this  example  are  carried  out  using  an 
IBM  1620  digital  computer.   Fig.  9  shows  the  dimensions  and  properties  of 
the  bridge  analysed  in  this  example.   Three  solutions  of  this  numerical 
example  were  accomplished.   The  first  and  second  solutions  used  the  consis- 
tent deformations  method  with  the  first  solution  taking  into  account  both 
the  effects  of  flexural  and  of  axial  stresses  and  the  second  solution 
considering  the  flexural  stresses  only.   The  third  solution  was  accomplished 
by  using  the  energy  method.  All  solutions  are  explained  separately  below. 


ARCH  EQUATION  :   y  =  _^_  (300  .  x). 


Ir   =  12IS  (ft4). 


A  =  I*  /3  (ft'), 
cr 


I*  =  rcr0o*ej  (Efc4)-   Aj  "  Tj  (ffc2>' 


E  =  30  x  10  psi. 


Fig.  9.  Dimensions  and  properties  of  a  girder-arch  bridge. 
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Solution  1. 

Considering  both  the  flexural  and  axial  stress  effects  and  using  the 
method  of  consistent  deformations. 

Solution: 

Following  the  steps  outlined  in  the  General  Analysis  Procedures  for 
solution  1,  the  final  results  are  obtained  and  shown  in  Fig.  10  to  Fig.  13, 
The  appropriate  data  are  shown  in  Tables  1  to  5  in  Appendix  A,   By  the  theory 
of  virtual  displacements,  it  can  be  seen  that  the  shape  of  each  influence 
line  is  reasonable  since  they  are  generally  of  the  shape  of  the  deflection 
curve.   This  solution  provides  the  most  accurate  results  since  it  takes  into 
account  the  most  stress  effects. 
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Fig. 10.  Influence  lines  for  the  redundant 

vertical  member  forces  in  solution  1. 
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11.  Moment  influence  lines  for  the  girder 
in  solution  1 . 
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Fig. 13 •  Final  influence  line  for  11   in  solution  1. 
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Solution  2. 

Considering  the  flexural  stress  effects  only  and  using  the  method  of 
consistent  deformations. 

Solution: 

Following  the  steps  outlined  in  the  General  Analysis  Procedures  for 
solution  2,  the  final  results  are  obtained  and  shown  in  Fig.  14  to  Fig.  17. 
The  appropriate  data  are  shown  in  Table  1  to  5  in  Appendix  B.   It  can  be  seen 
that  most  of  the  results  are  quite  close  to  those  of  solution  1.   It  is  pre- 
dictable that  the  results  for  the  influence  line  for  X.(or  R. )  would  not  be 
very  accurate  since  the  actual  structural  behaviour  for  this  member,  which  is 
the  longest  among  the  vertical  truss  members,  is  quite  far  removed  from  the 
assumption,  of  negligible  axial  stress.   It  may  be  noted  that  because  of  the 
assumptions  made  and  the  node  points  selected  there  is  no  information  on  the 
moment  influence  lines  for  the  girder  and  arch  for  the  segments  between  node 
points  A,  1  and  11,  B. 
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Kig.l'l.  Influence  lines  for  the  redundant  vertical 
member  forces  in  solution  2. 
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Fig. 15«  Moment  influence  lines  for  the  girder 
in  solution  2 . 
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Fig. If).  Moment  influence  lines  for  the  arch 
in  so iu t ion  2 . 
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Fig. 10.  (Contd. ) 


Fig. 17 •  Influence  line  for  H 


in  Solution  2 , 
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Solution  3.  Energy  method  solution. 

Solution: 

Following  the  steps  outlined  in  the  General  Analysis  Procedures  for 
solution  3,  the  final  results  are  obtained  and  shown  in  Fig.  18,  Fig.  19  and 
Fig.  20.   The  appropriate  data  are  shown  in  Table  1  and  Table  2  in  Appendix 
C.   It  may  be  seen  that  the  results  are  quite  close  to  those  of  solutions  1 
and  2.   Since  in  this  method  of  solution,  the  vertical  member  forces  are  not 
relevant  there  is  no  information  obtained  concerning  them.   Checking  the  re- 
sults, it  can  be  seen  that  the  moments  at  points  1  and  11  are  quite  far  re- 
moved from  the  values  obtained  from  solutions  1  and  2.   This  can  be  explained 
by  the  fact  that  the  actual  structural  behaviour  differs  a  great  deal  from 
the  assumed  behavior.   Points  1  and  11  are  hinged  for  the  arch  while  the  cor- 
responding points  in  the  girder  are  continuous.  Therefore  the  assumption 
that  the  curvature  due  to  load  in  the  arch  is  equal  to  the  curvature  due  to 
load  in  the  girder  is  in  reality  a  relatively  poor  assumption. 
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inal  influence  line  for  H   in  solution  3. 
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A  comparison  between  the  results  of  solution  1,  2,  and  3  shows  that 
they  are  reasonably  close  to  each  other.   The  first  solution  provides  the 
most  accurate  results  but  requires  the  most  computational  effort.   The  second 
solution,  which  requires  less  computational  effort,  provides  very  accurate 
results  in  comparison  with  those  of  solution  1.   The  last  solution  provides 
the  most  inaccurate  results.  They  are,  however,  still  reasonably  close  and 
require  considerably  less  computation  than  either  of  the  other  solutions. 
Of  these  three  solutions  the  third  solution  might  well  be  used  for  a  prelimi- 
nary design  while  either  of  the  other  two  might  be  used  for  the  final  design. 
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CONCLUSIONS 

By  comparing  the  three  sets  of  solutions  of  the  numerical  examples,  it 
can  be  seen  that  the  results  correspond  closely. 

For  solution  1  the  concept  of  consistent  deformations  which  is  very 
familiar  to  structural  engineers  was  used.  With  the  consideration  of  thrust 
in  the  arch  and  the  axial  stresses  in  the  vertical  members,  the  procedures  of 
computation  are  much  more  troublesome  than  they  are  in  the  other  two  solutions. 
However,  this  is  the  most  accurate  solution  of  the  three.  Solution  2  makes 
use  of  the  same  concept  of  consistent  deformations  but  includes  the  assumption 
that  the  thrust  and  axial  stress  effects  can  be  ignored.   The  results  are 
nearly  the  same  as  those  of  the  first  solution.   Hence  it  would  seem  that  It 
is  satisfactory  to  use  this  solution  in  a  practical  design.   In  the  last 
solution,  which  makes  use  of  the  general  concept  of  minimum  energy  and  the 
theory  of  virtual  displacements,  the  redundants  are  greatly  reduced  from  14 
degrees  to  3  degrees.   This  reduction  of  redundants  is  a  great  advantage  in 
the  computations  and  makes  it  practical  to  carry  them  out  by  hand  computations 
in  case  a  computer  is  unavilable.  As  a  matter  of  fact,  the  energy  method  can 
be  used  in  the  preliminary  design  without  resulting  in  significant  error. 

It  should  be  noted  that  the  moment  influence  lines  for  point  1  and  2  for 
the  energy  method  are  quite  far  removed  from  those  obtained  from  solution  1 
and  solution  2  by  the  consistent  deformations  method.   This  discrepancy  can  be 
explained  by  the  assumption  made  in  the  energy  method  solution  that  the  angle 
changes  along  the  girder  and  the  arch  rib  are  equal.   It  can  be  seen  that 
actually  the  arch  Is  hinged  at  point  c  which  corresponds  to  point  1  of  the 
girder  which  is  continuous  at  that  point.   That  is,  the  assumption  of  equal 
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change  of  angles  along  the  girder  and  the  arch  at  the  corresponding  points  of 
the  structure  deviates  significantly  from  the  actual  structural  behaviour. 
It  is  not  therefore  surprising  that  the  various  solutions  of  the  problem  also 
have  significant  deviations  at  these  points.   If  the  arch  portion  were  a 
fixed-ended  arch,  the  differences  between  the  results  should  be  greatly 
decreased. 
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NOTATION 
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4 
is 

ir 
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m , ,  m,  -  etc. 


mij'  mij 


Mtj'  "ij 


'kt 


Pj»  Pk»  etc- 


*U 


Cross  sectional  area  of  arch  rib. 

Cross  sectional  area  of  vertical  hinged  members. 

Modulus  of  elasticity. 

Horizontal  reaction. 

Final  horizontal  reaction  of  the  girder  arch  bridge  due  to 
loads  acting  at  node  point  j. 

Horizontal  reaction  of  a  single  two-hinged  arch  due  to  loads 
acting  at  node  point  j  of  a  two-hinged  arch. 

Moment  of  inertia  at  the  crown  of  arch. 

Moment  of  inertia  of  girder. 

Moment  of  inertia  of  arch. 

Length  measured  in  the  X-axis  direction. 

Length  of  the  vertical  member  j. 

Moment  caused  by  Vc,  VD,  and  H  in  the  reduced  structure  In 

the  energy  method  respectively. 

Moment  due  to  loads  acting  at  node  point  j,  k,  etc. 

Moment  induced  at  node  point  1  caused  by  a  unit  load  acting 
at  node  point  j  for  girder  and  arch  respectively. 

Moment . 

Total  induced  moment  in  girder  and  arch  respectively. 

Total  induced  moment  at  section  k  of  the  compound  section 
of  girder  and  arch  in  the  energy  method. 

Thrust  due  to  loads  acting  at  node  points  j,  k,  etc. 

Loads  acting  at  node  points  j,  k,  etc. 

Reaction  at  point  1  due  to  loads  acting  at  point  J. 

Length  measured  along  arch. 
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U  :   Strain  energy. 

vc»  VQ        :   Vertical  reactions  at  hinged  point  C,  D  respectively. 

X,,  X  ,  etc.   :   Forces  in  the  vertical  truss  members  i,  j,  etc. 

Xj,  :   Forces  induced  in  vertical  member  i  due  to  loads  acting  at 

node  point  j . 

X  :   Length  of  segments  between  adjacent  node  points  measured  in 

the  X-axis  direction. 

Vertical  deflections  of  girder  at  section  k,  j ,  etc. 

Shortening  or  elongation  of  the  vertical  members  k,  j,  etc. 

Vertical  deflection  of  arch  at  section  k,  j ,  etc. 

Deflection  at  section  k  due  to  loading  acting  at  node  point  j. 

Deflection  at  section  k  of  a  simple  girder  with  loads  acting 
at  node  point  j. 

Deflection  at  section  k  of  a  simply  supported  arch  with  loads 
acting  at  node  point  j  of  the  simply  supported  arch. 

Vertical  displacement  of  hinged  points  C  and  D  respectively 
in  the  energy  method. 

Relative  horizontal  displacement  of  hinged  points  C  and  D  In 
the  energy  method. 

Angle  between  the  X-axis  and  the  tangent  of  the  arch  at  point 
j. 
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ABSTRACT 

The  consistent  deformation  method  and  the  energy  method  are  used  In 
this  report  for  the  stress  analysis  of  a  girder-arch  bridge.   Using  the 
first  method,  two  solutions  were  calculated,  the  first  of  which  takes  into 
account  the  axial  stress  effects  in  the  structure,  and  the  second  ignores 
these  effects.  A  third  solution  was  developed  using  the  energy  method. 

Of  the  many  possible  methods  of  analysis  for  this  type  of  bridge,  this 
report  presented  two  which  may  be  used  to  analyse  high-order  statically 
Indeterminate  structures.  A  digital  computer  was  available  and  therefore 
the  calculations  for  both  of  these  methods  were  wholly  performed  by  the 
digital  computer.  Since  the  bridge  analysed  in  this  report  has  variable 
cross  sections  in  the  arch,  this  report  shows  a  general  procedure  for  a 
computer  analysis  to  handle  a  structure  with  variable  cross  sections. 

Numerical  examples  are  given  to  illustrate  the  procedures  for  finding 
the  Influence  lines  for  the  structure.   Since  three  solutions  were  accom- 
plished for  the  same  example,  a  comparison  of  the  results  of  the  various 
methods  can  be  made.  The  comparison  of  the  three  solutions  Indicates  that 
the  results  are  reasonably  close  to  each  other. 


